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We present a general, analytic recipe to compute the entanglement that is generated between arbitrary, discrete 
modes of bosonic quantum fields by Bogoliubov transformations. Our setup allows the complete characterization 
of the quantum correlations in all Gaussian field states. Additionally, it holds for all Bogoliubov transformations. 
These arc commonly applied in quantum optics for the description of squeezing operations, relate the mode 
decompositions of observers in different regions of curved spacetimes, and describe observers moving along non- 
stationary trajectories. We focus on a quantum optical example in a cavity quantum electrodynamics setting: 
an uncharged scalar field within a cavity provides a model for an optical resonator, in which entanglement is 
created by non-uniform acceleration. We show that the amount of generated entanglement can be magnified 
by initial single-mode squeezing, for which we provide an explicit formula. Applications to quantum fields in 
curved spacetimes, such as an expanding universe, are discussed. 
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1. Introduction 

Over the past decade the discipline of relativistic quantum information has received much at- 
tention (see Ref. [1-7] for a selection and Refs. [8, 9] for reviews). Its aim is the study of the 
resources and tasks of quantum information science in the context of relativity. In particular, 
finding suitable ways to store and process information is a main goal. 

It has therefore been a central focus of previous efforts to identify the degradation effects due 
to relativistic, accelerated motion [10-16] and spacetime curvature [17, 18]. At the centre of 
these deteriorations lie two ingredients: the Bogoliubov transformation of the mode operators 
and the presence of a horizon. While the Bogoliubov transformations are responsible for the 
creation and shift of excitations, the horizon causes some of these excitations to be inaccessible 
to the observer. This usually leads to a loss of information and, consequently, a degradation of 
the quantum correlations. A rare exception of a situation where entanglement is generated by 
acceleration can be found in Ref. [19]. 

In Ref. [20] cavities were proposed as a suitable way of storing and processing information in the 
relativistic context. Cavities, represented by appropriate boundary conditions, can be uniformly 
or non-uniformly accelerated without creating a horizon. Nevertheless, the absence of an event 
horizon does not guarantee complete access to all modes of a quantum field. In particular, it 
was shown in Refs. [21, 22] how the entanglement between a non-uniformly moving cavity and 
an inertial reference cavity is degraded when only particular modes are being considered. In 
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this case the information loss can be prevented by appropriately timing the periods of uniform 
acceleration and inertial motion. 

However, the role of the Bogoliubov transformations lies foremost in the generation of en- 
tanglement, an effect that can be found for a wide variety of situations. In quantum optics 
they are used for the description of single- mode as well as multi-mode squeezing operations, 
see, e.g., Ref. [23]. Entanglement is also created between modes of quantum fields in curved 
spacetimes [24, 25]. And, recently, entanglement generation has been demonstrated and studied 
quantitatively for non-uniform cavity motion in flat Minkowski spacetime in Refs. [26, 27]. 

Here we present a framework that generalises and significantly simplifies the previous ap- 
proaches on entanglement generation in relativistic quantum information. We employ continuous 
variable techniques [28] that allow us to investigate the quantum correlations that are created by 
arbitrary Bogoliubov transformations of any discrete number of modes. The transformed state 
of these modes is described by its covariance matrix, which completely determines the entan- 
glement of the system. Results for any given Bogoliubov transformation can thus be obtained 
analytically in principle, if the transformation coefficient are known. 

We apply our method to a quantum optical setting: an optical cavity that is moving non- 
uniformly is modelled by confining an uncharged scalar field by appropriate boundary condi- 
tions [21]. The cavity is following a trajectory that consists of segments of uniform acceleration 
and inertial motion. We work in a perturbative regime where the products of the cavity's width 
and the individual accelerations in all segments are small, but all calculations can be carried out 
analytically. The Bogoliubov coefficients can be expressed as a Maclaurin series in h, a parame- 
ter which represents the product of the acceleration at the centre of the cavity and the cavity's 
width S. The creation of entanglement for initially uncorrelated, single-mode squeezed states of 
the field is quantified in terms of the expansion coefficients of the power series. We find that the 
entanglement generation in this fully relativistic setting can be enhanced by initial single-mode 
squeezing. This provides an additional controllable parameter for possible experimental setups 
to enhance the visibility of the effect. 

We further emphasize the applicability of our approach beyond cavity quantum electrodynam- 
ics by verifying the results obtained for an expanding universe in Ref. [24] as a curved spacetime 
examples. 

This article is structured as follows. In Sec. 2 we establish the basic description of Bogoliubov 
transformations for arbitrary Gaussian initial states in terms of the corresponding symplectic 
transformation of the covariance matrix. In Sec. 3 we revise the description of non-uniformly 
moving cavities from Refs. [21, 22] and, subsequently, apply the formalism of Sec. 2 to study the 
entanglement creation by motion in this scenario. 

2. Bogoliubov transformation of Gaussian states 

It is the aim of this section to present a general machinery that allows a simple characterisation of 
the mode entanglement of bosonic quantum fields. In particular we aim to construct a framework 
that works for any specified Bogoliubov transformation. 

To this end, let us consider an arbitrary discrete set of mode functions {4> n \ n = 1, 2, 3, . . .} of a 
bosonic quantum field, let us assume a scalar field for simplicity, with associated annihilation and 
creation operators, a n and a) n , respectively. The field operators satisfy the canonical commutation 
relations [a m ,a n ] = [aln^al] = and [a m ,On] = 8 mn . The functions <j) n , which we take to be the 
solutions to a (relativistic) field equation, e.g., the Klein-Gordon equation, form a complete set 
of orthonormal modes with respect to a chosen scalar product, i.e., (4> m ,(p n ) = 5 mn , see, e.g., 
Ref. [29]. In the usual Fock representation the vacuum state is annihilated by all a n , i.e., a n | ) = 
0, Vn, while particle states are created by the (repeated) action of the creation operators a$i- 

We can then perform a Bogoliubov transformation that relates our initial modes <fi n 
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and the associated operators a n and a n to another (complete, orthonormal) set of modes 
{4> n | n = 1, 2, 3, . . .}, with operators a n and a„, respectively, i.e., 



Here we have used the conventions of Ref. [29]. The asterisk denotes complex conjugation and 
«mn = (f^m, <^n) and (3 mn = — (</> n , 4>* m ) are the coefficients of the Bogoliubov transformation. 

This transformation can represent a physical operation in a laboratory, e.g., a squeezing op- 
eration in an optical setup, or a change of potential in a harmonic oscillator chain of trapped 
atoms. Alternatively, the transformation may describe a change of observer in a (curved) space- 
time, e.g., a black hole spacetime or a cosmological model, or relate the modes at different times 
of a single observer that is moving non-uniformly. For a continuous spectrum of solutions the 
sum in Eq. (1) has to be replaced by an appropriate integral, but we will specialise to discrete 
sets of modes in the following. 

Analogously to before, the vacuum state | ) of the new mode decomposition is defined by the 
property a n | ) = 0, \/ n. However, an important feature of the Bogoliubov transformations is the 
non-equivalence of the corresponding vacua, i.e., a m | ) =0 only if j3 mn = 0, Vn. Consequently, 
the vacua generally have a nontrivial relationship, see, e.g., Ref. [30]. 

To study the quantum correlations between particular modes two more steps are required. 
The complete Fock state basis of the original modes has to be expressed in terms of the Fock 
states of the transformed modes. Subsequently, the subset of modes that is not considered must 
be traced over. While the first step can generally be accomplished, the tracing procedures are 
usually cumbersome. 

In the following we demonstrate how these computations can be simplified significantly. Let 
us switch to the covariance matrix formalism discussed in Ref. [28]. This means, instead of 
the complete Fock space description of the quantum state, we consider only the matrix a with 
components 

aij = ( XiXj + XjXi ) - 2 ( X< ) ( Xj ) . (2) 

The quadrature operators Xj are chosen to be the generalised positions and momenta, i.e., 
X(2n-i) = ^jffan + a «) an d X(2n) = ^|(on — , where the index n = 1,2,3, .. . labels the 
modes, and ( O ) denotes the expectation value of the operator O. The real, symmetric covariance 
matrix a together with the vector of first moments Xj completely characterises all Gaussian 
states, i.e., states that can be described by quasi-probability distributions of Gaussian shape in 
phase space, see, e.g., Ref. [28]. However, the covariance matrix itself is a sufficient description 
of all properties pertaining to entanglement. While these continuous variable tools have been 
frequently employed in quantum optics settings, their application to quantum field theory and 
relativistic quantum information was previously explored only in limited situations [12, 17]. 

Unitary transformations in the Fock space are represented by symplectic transformations in 
phase space. A transformation S is called symplectic, if it leaves the symplectic form invariant, 
i.e., SQS T = $7, where [Xi, Xj] = and {Sij) T = (Sji). From Eq. (1) it is straightforward to 
find the expression for the symplectic representation S of a Bogoliubov transformation in terms 
of its general coefficients a mn and j3 mn . The matrix S is decomposed into 2x2 blocks A4 mn as 



(Mu M12 M Vi . . \ 

M 2 i M22 M23 ■ ■ 

M u M32 M33 ■ ■ 

V i i ; 



(3) 
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where the sub-blocks are given by 

Here 3ft(z) and ^(z) denote the real part and imaginary part of the complex number z respec- 
tively. The transformed covariance matrix a is then simply obtained as 

a = SaS T . (5) 

The crucial technical simplification lies in the ensuing step. The partial trace over any subset 
of modes £ C {n = 1, 2, 3, . . .} can be computed trivially by eliminating all rows and columns 
corresponding to modes m £ £ . This technique applies for any initial state a. However, if 
we restrict ourselves to Gaussian states, the unitarity of the transformation ensures that a 
remains Gaussian, which allows us to use the specifically useful tools for the quantification of 
Gaussian entanglement [31]. In particular, we are interested in the generation of entanglement 
from initially uncorrelated states. In this case the covariance matrix a has a block diagonal 
structure, a = diag{^ n }, where the 2x2 blocks ij) n are locally equivalent to single-mode squeezed 
states with squeezing parameters s n , i.e., 

^ = (v A) • (6) 

Let us now consider the transformed state of two modes, labelled k and k' respectively. In other 
words, £ = {n = 1, 2, 3, . . . \n ^ k, k 1 } is the set of all modes except the chosen pair k and k' . 
The transformed covariance matrix dkk' of these two modes is then given by 

/ Cfcfc Ckk'\ 

°kk' = [ „ n ( 7 ) 



Ck'k Ck' 



where dj = £ n M in ip n Mj n . 

Let us now turn to the quantification of the entanglement that is generated between the modes 
k and k' . The quantity that encodes all the information about the entanglement of symmetric 
two-mode Gaussian states is the smallest symplectic eigenvalue P_ of akk' — ?V &kk' , where 
Tfc/ = diagjl, 1, 1, — 1} represents the partial transposition of mode k' , see Ref. [28]. For this 
class of states the Gaussian measures of entanglement, such as the (logarithmic) negativity or 
the Gaussian entanglement of formation, are monotonously decreasing functions of z?_ . However, 
the two-mode Gaussian state dkk' of Eq. (7) is generally not symmetric with respect to the 
interchange of k and k' . Consequently, different entanglement measures can exhibit different 
ordering of two entangled states. 

With this in mind we specialise to the use of the negativity M as an entanglement measure. 
It has the advantage of being easily computable, i.e., the negativity of a^k' is given in terms of 
the smallest symplectic eigenvalue P_ of dkk' by the simple formula 

AA = max{0, (1 -VJ)/2u-}. (8) 

Furthermore, it can be easily compared to results obtained for non-Gaussian states in Ref. [26]. 
The smallest symplectic eigenvalue is obtained by diagonalising the matrix Okk' by a symplectic 
operation D, such that Dokk 1 D T = diag{— z?_, — z?+,z?+} , where < P_ < p + . The quanti- 
ties ±P± can be computed as the eigenvalues of i^ldkk' i n a straightforward way. For < V- < 1 
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the state &kk' is entangled. 

Other well known entanglement measures for Gaussian states, e.g., Gaussian entanglement 
of formation, can also be calculated straightforwardly for two-mode states in principle [28]. 
The transformed state will generally be a mixed state but all calculations can be done analyt- 
ically, provided that the Bogoliubov coefficients are given and the infinite sums in Eq. (7) are 
convergent. In Sec. 3 we discuss an example of a discrete spectrum for which the Bogoliubov co- 
efficients are given as a perturbative expansion and the sums converge [21, 22]. We compute the 
entanglement generated from non-uniform cavity motion for particular initial states of interest. 
For continuous spectra, however, the corresponding integrals are often known to be divergent, 
see, e.g., Ref. [32], thus rendering any conclusive statements about the entanglement generated 
between the modes fruitless, unless the coefficients have a simple structure. 

Such a situation presents itself in the case of a charged, scalar field in an expanding universe, 
which was discussed in Ref. [24]. There the spectrum of the quantum field is continuous, but 
the Bogoliubov transformation between the asymptotically flat remote past and future cou- 
ple only modes of opposite momenta. In this fashion, the transformation provides an effective 
discretization and we can reproduce the results of Ref. [24] with our methods. 



3. Non-uniform cavity motion 

An example for a discrete, bosonic spectrum is obtained by confining a scalar quantum field 
to a cavity by appropriate Dirichlet boundary conditions in (1 + 1) dimensions. As proposed 
in Ref. [21] the rigid cavity can follow a worldline that is composed of segments of inertial 
motion and uniform acceleration. This non-uniform motion generates non-trivial Bogoliubov 
coefficients, which result in a generation of entanglement between the modes inside one cavity [26] 
and, consequently, lead to a degradation of initial entanglement between modes in different 
cavities [21]. 

For a quantitative description the Bogoliubov transformations are expanded as a Maclaurin 
series in the small, dimensionless parameter h, i.e., 

a = + a™ + a™ + 0(h 3 ) , (9a) 
P = /?« + 0(2) + 0(h s ) , (9b) 

where the superscripts ^ denote quantities proportional to h n . The parameter h is the product of 
the cavity's length in its instantaneous rest frame and the proper acceleration at the centre of the 
cavity. Here we use units such that Planck's constant and the speed of light are dimensionless 
constants, h = c = 1, and the entity 0(x)/x is bounded as x goes to 0. Furthermore, the 
coefficient must include the phases of the free time evolution in the uniformly accelerated and 
inertial segments, while it reduces to the identity for vanishing accelerations. We therefore have 
ctmn = 5 rnn G m , where G m is a mode-dependent phase factor of unit magnitude, i.e., \G m \ = 1. 
Additionally, we find that the linear corrections vanish on the diagonal, i.e., cxnn = Pnn = 0. 

The perturbative calculations require a closer inspection of the techniques for the calculation of 
the symplectic eigenvalues P_ from Sec. 2. Let us assume that the initial state Ukk 1 is transformed 
by the Bogoliubov transformation according to 

&kk' = &kk' + <y C kk' > (1°) 

whereof, is a small correction to the initial state. For uncorrelated, pure initial states <7fcfc'i 
which we want to study here, the unperturbed symplectic eigenvalues of T# are 
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degenerate, i.e., u\ = v_ = 1. This requires the diagonalisation of the subspaces of the degen- 
erate eigenvalues of the perturbation a kk , to obtain the corrected smallest symplectic eigenvalue 

P_ = 1 ± z/f. In other words, the corrections ±z/f to = 1 are given by the eigenvalues of 
the matrix (7^) with components 

lij = (e i± \in a c kk , I e j± ) , (i, j = 1,2), (11) 

where \&j±) is the j-th eigenvector of i£lo kk > corresponding to eigenvalue ±1. An important 
feature of the unperturbed state o~ kk >, an d consequently of the eigenvectors of i$la kk > > is its 
time dependence due to the free time evolution of the modes, i.e., o~ kk > ~~ ^ Ro~kk'R T ■ The local 
rotation R is represented by an orthogonal, block diagonal matrix R = diag{i? n }, where the 
block of the n-th mode can be written as R n = 5i (G n )t2 + i ^(G n )ay, and a y is the usual Pauli 
matrix. The leading order correction to the negativity (8) is given by |z^£|/2. As a particular 
example we can study the case of symmetric, initial single-mode squeezing in the modes k and 
k! , with (real) squeezing parameters s k = s k < = s (see Eq. (6)). In this situation the leading 
order corrections to the negativity become 

M = (R(Gft$) a + ( HGl^l) cosh( s ) - 9f(GJagJ) sinh( S )) 2 ) 1/2 + 0(h 2 ) . (12) 

For modes of opposite parity, i.e., if (k + k') is odd, the linear coefficients are non-zero see 
Ref. [21]. In the case of vanishing squeezing parameters, s k = s k < = 0, the expression in (12) 

reduces to Af = \/3 k 2\, which is consistent with the expression for the entanglement that is 
generated from the bosonic vacuum in the Fock representation, see Ref. [26]. As can easily be 
seen from Eq. (12) the initial single mode squeezing introduces a non-negative term into the 
negativity. Thus the generated negativity is always enhanced with respect to the vacuum case. 
An example for the effect of single-mode squeezing on the generation of entanglement is shown 
in Fig. 1(a) for a particular travel scenario of the cavity. It can be readily observed, that the 
negativity grows with e s for e s 3> e~ s . While this fact can be utilised to enhance the visibility of 
the entanglement generation, it also limits the validity of the perturbative regime, which requires 
further investigation. 

A transparent method to quantify the perturbation of the initial state is the analysis of the 
system's mixedness. As discussed in Ref. [21, 27] the Bogoliubov coefficients that relate the 
modes k and k' can be consistently renormalized to represent a unitary transformation on the 
subspace of these two modes. This two-mode truncation, which leaves the linear coefficients ctili 

and pj^, invariant, can be implemented if k and k' have opposite parity. The state of a truncated 
system of this type will undergo a unitary evolution due to the cavity motion. In particular, 
the determinants of the individual mode subspaces, det Ckk and det Ck'k' , will be identical to 
leading order in h, i.e., the state is symmetric. 

We then exploit the fact that any pure, symmetric two-mode Gaussian state is locally equiva- 
lent to a two-mode squeezed state [28]. This enables us to express the corresponding two-mode 
squeezing parameter r as a Maclaurin series in h by using |r| = ^ arsinh y/— det C^k 1 - We recover 
the expression of Eq. (12) for the negativity, i.e., M = |r| as is expected for two-mode squeezed 
states. This suggests that Eq. (12) is valid as long as the perturbation to the mixedness of 
the system is small, i.e., as long as |detcjfcfc' — 1| <C 1. We compute the determinant of the 
transformed state for the symmetrically single-mode squeezed initial state to be 

det d kk , = 1 + 4(it fc , + / fe %)(cosh S + 1) + 4(/ fe %, + f^ k )(coshs - 1) 

-4sinh, £ ^a^> + a^>), (13) 

n^k,k' 
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Figure 1. The leading order contribution Af/h to the negativity, see Eq. (12), is shown in Fig. 1(a) for a massless scalar 
field confined to a moving cavity of width 8. The cavity's trajectory is inertial outside a single segment of uniform, linear 
acceleration h/5 of duration r, as experienced at the centre of the cavity. The Bogoliubov coefficients for this travel scenario 
can be found in Rcf. [21]. Plots are shown for an initial state of the chosen modes (k,k') = (1,2) that is symmetrically 
single-mode squeezed with s = 1 (dashed, red) or in the (displaced) vacuum, s = (solid, blue). Fig. 1(b) illustrates the 
behaviour of the quantities F k k i /h 2 = (/^ fe / + fy^^/h 2 ^ or i ncreasm g mode numbers. The same field and travel scenario 
as in Fig. 1(a) are used, but the mode numbers are varied. From bottom to top the curves show (k,k') = (1,2) (solid, 
green), (k,k') = (10,11) (dashed, red), (k,k') = (1,10) (dashed, blue), (fc, k') = (20,21) (solid, red) and (fc, k') = (1,20) 
(solid, blue). All curves are plotted as functions of the temporal parameter u := hr/[AS atanh(h/2)]. 

where f^ k , = ^En^<l«2l 2 and fLv = 3 En^fc' itfftl 2 - The "-coefficients grow with the 
mode numbers of the selected modes [21], while \(3mn\ — > if m — > oo or n — > oo. We can thus 
constrain the range of validity of the perturbative result of Eq. (12) by the simple requirement 
-Ffc,fc'(coshs-l) < 1, where F k>k < = (f k ^ k ,+f k ,^ k )- Sample plots for F k ^/h are shown in Fig. 1(b). 
Additionally, when a massive, bosonic field is considered the mass of the field excitations must 
be appropriately restricted, see Ref. [21]. However, as long as these conditions are satisfied, the 
single- mode squeezing parameter, the choice of modes and the particular travel scenario for the 
cavity can be tuned to enhance the entanglement generation effect. 

4. Conclusions 

We have presented a general framework that allows the quantification of the entanglement that 
is generated between arbitrary modes of bosonic quantum fields by Bogoliubov transformations. 
Our setup combines techniques from quantum optics, quantum field theory and quantum infor- 
mation procedures to describe the quantum correlations that arise from non-uniform motion, 
spacetime curvature and quantum optical operations in a covariance matrix formalism. This 
removes the necessity for cumbersome partial tracing procedures to quantify the entanglement 
between field modes. For Gaussian initial states we can fully characterise the entanglement that 
is produced in these situations. 

We have discussed a particular example, a cavity containing a relativistic quantum field, 
in which non-uniform motion was recently found to create entanglement [26]. The Bogoliubov 
coefficients for this scenario are given as a perturbative expansion and we compute the negativity, 
the entanglement measure of our choice, to leading order in the expansion parameter. We discuss 
in detail the regimes in which the perturbative calculations can be trusted. 

To leading order in the expansion parameter, the transformed state of two modes of opposite 
parity is equivalent to a pure, two-mode squeezed state. The corresponding squeezing parameter 
can be used as an alternative route to quantify the entanglement of the system. We give clear 
criteria for the validity of this two-mode truncation of the Bogoliubov transformation. The 
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calculation of the mixedness of the two-mode state without the truncation provides a simple 
benchmark against which to judge the limits of the small parameter expansion. 

By application of the methods presented here we are able to extend the previously found 
results to the important class of Gaussian states. We find that initial single- mode squeezing can 
enhance the correlations that can be generated from the vacuum, thus greatly improving the 
prospects of experimental verification of this effect and of its implementation in new quantum 
technologies. 

Under some circumstances, the approach of this article even permits application to particular 
cases of continuous spectra of quantum fields, such as a charged, scalar field in an expanding 
universe [24]. 
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